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Abstract. We study a deformation of the nonlinear Schrodinger equation recently 
derived in the context of deformation of hierarchies of integrable systems. This 
systematic method also led to known integrable equations such as the Camassa-Holm 
equation. Although this new equation has not been shown to be completely integrable, 
its solitary wave solutions exhibit typical soliton behaviour, including near elastic 
collisions. We will first focus on standing wave solutions, which can be smooth or 
peaked, then, with the help of numerical simulations, we will study solitary waves, 
their interactions and finally rogue waves in the modulational instability regime. 
Interestingly the structure of the solution during the collision of solitary waves or 
during the rogue wave events are sharper and have larger amplitudes than in the 
classical NLS equation. 
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1. Introduction 

In this work, we will study a deformation of the nonlinear Schrodinger equation, called 
hereafter the CH-NLS equation. The standard dimensionless form of this equation is 
given for a complex field u by 

imt + Uxx +{\u\^ — a^\ux\^) = m = u — a‘^Uxx, (1) 

where cr = ±1 selects between the focussing and defocusing equation. Although the two 
equations have different properties, we will only consider the focussing case, thus we 
will set a = 1 for the rest of this paper. 

This equation was recently derived by the author in [T] when developing a theory of a 
deformation of hierarchies of integrable systems. When applied to the AKNS hierarchy 
of |2] this deformation led to the derivation of several known integrable equations. 
The most important one is the Camassa-Holm equation [3], found as a deformation 
of the KdV equation. The CH equation has been extensively studied in the context 
of shallow water waves, see 011]. This has led to the development of the a-model, a 
regularized fluid model which extends the CH equation to higher dimensions (see [5] and 
references therein). Other deformations of integrable equations derived in [1] include the 
modified CH equation j6| and other coupled equations [7| viewed as deformations of the 
mKdV and coupled KdV equations respectively. These equations all correspond to the 
deformation of the KdV flow of the AKNS hierarchy. However, the method of [T| can 
deform the previous flow as well, which corresponds to the NTS equation. The CH-NLS 
Equation ([^ is the result of this systematic approach which gathers several integrable 
equations together in a single framework. Although most of the equations derived in this 
context were independently shown to be completely integrable, this deformation does 
not guarantee their complete integrability. For this reason, there is still no proof for 
the complete integrability of the CH-NLS Equation ([^. The question of the complete 
integrability of the CH-NLS equation therefore remains an open problem. 

The main focus of this paper we will be on the solutions of this equation, which 
share the same properties as the solutions of integrable equations. We found a one 
parameter family of standing waves solutions with the remarkable property of having 
a transition between smooth and peaked shapes by varying this parameter. Although 
very rare in dispersive equations, such peaked solutions have already been found and 
studied in a generalised NLS equation by 0|9|. We then numerically studied solitary 
waves by running initial value problems with hyperbolic secant initial condition^ We 
observed the emergence of dispersive and solitary waves. The latter can even have 
a breathing behaviour, typical of cubic nonlinear equations such as mKdV or NLS for 
example. Numerical simulations of the collisions between these solitary waves show that 
they keep their identity after a collision, despite numerical errors. The existence and 
the stability of solitary waves in non-integrable equations is a vast and difficult subject 
which will not be addressed here. We refer to [lO] for a recent review. 

f Videos of some simulations are available on http://wwwf.imperial.ac.uk/-aal0213 
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It is well known that the existence of solitary waves in most of the nonlinear partial 
differential eqnation relies on the balance between nonlinearities and linear dispersions. 
One of the few exceptions is the Camassa-Holm eqnation [3] which admits solitons 
even in the dispersionless case. The fnndamental mechanism of the CH eqnation is the 
nonlinear dispersion, coming from the Helmholtz operator. By construction, the CH- 
NLS Equation Q has a similar nonlinear dispersion based on the Helmholtz operator 
but no transformation can remove its linear dispersion, as it is usually used for the study 
of the CH equation. The nonlinear dispersion based on the Helmholtz operator is in fact 
a common feature of CH type equations, a growing family of integrable equations which 
admit soliton solutions without linear dispersion (see [sum El Cl [HE] among others). 
The CH-NLS equation is thus a very particular equation where the linear dispersion 
plays a more important role than for other equations of this type. 

Another interesting fact is that if m is considered as the fundamental held then u is a 
nonlocal function of m through the relation u = K * m with K = The linear 


dispersion relation is of the form u{k) = 


which is a bounded function. Waves of 


high frequencies with respect to a will thus not propagate. This effect can be observed 
in the simulations of collisions where high frequency waves were artihcially created by 
the numerical scheme and then stayed at the same position. This type of modihcations 
of the linear dispersion relation for the NTS equation has already been mathematically 
investigated for example by na in the context of short pulses in optical hbres. We 
will also see the CH-NLS equation shows more extreme behaviours than the standard 
NLS equation, in particular during collisions of solitary waves or rogue waves in the 
modulational instability regime. 

In Section we will look at the conservation laws of the CH-NLS equation as 
well as a short discussion on its non integrability. The peaked and smooth standing 
wave solutions will be presented in Section the numerical initial value problems 
meant to isolate and characterise solitary waves are in Section and their collision 
in Section jS] The collisions exhibit more extreme behaviours than the collisions of NLS 
solitons, namely a higher amplitude and sharper structures during the collisions. The 
same features are found in the modulational instability, studied in the last Section 
where the CH-NLS Peregrine type solutions are higher and sharper than the NLS 
Peregrine solutions. Appendix A contains the numerical scheme used in this work and 
a convergence analysis for one of the collision of solitary waves. 


2. Conservations laws and symmetries 

The CH-NLS equation is a Hamiltonian system with Hamiltonian given by 


P = - 
2 


{mUx — mUx)dx 


and with the non-canonical Hamiltonian structure of the form 

—2md~^m dx + 2md~^fn 
dx + 2fnd~^m —2fnd~^m 


K = 


( 2 ) 


(3) 
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This Hamiltonian structure is exactly the same as second NLS Hamiltonian structure. 
The Hamiltonian P corresponds to the momentum of the solution, as it is associated to 
the space translation invariance of the equation. As will will discuss at the end of this 
section, the Hamiltonian corresponding to the time invariance of this equation has not 
yet been found, if it exists. A similar feature already exists as noted in [Tn| where the 
same conserved quantity arises form the space and the time invariance of the equation 
independently. 

A essential symmetry of the NLS equations is their invariance under phase shifts, 
which is also shared by the CH-NLS equation. For the NLS equation this symmetry is 
associated with the conservation of the total mass 

M= \m\‘^dx (4) 


which is also a conserved quantity of the CH-NLS equation. The associated flow is given 
by = Jq (|^, = {im, —ifnf'^ where Jq := J{a = 0) from Equation (6). 

However, as we will see below, this Hamiltonian structure does not produce the CH-NLS 
equation and the link between the conservation of mass and the phase shift symmetry 
seems broken, despite the existence of both properties. 

The CH-NLS equation is invariant under the following “scaling transformation” 

f —)■ a; —)■ u{t,x) ^ Xu{X‘^t, Xx) and a —)■ A“^a. 


This is not a proper transformation as it only relates solutions with different values of 
a. Nevertheless this scaling can be useful to obtain new solutions and makes clear the 
importance of a to hx a particular length scale in this equation. 

Another important symmetry of the NLS equation is the Galilean invariance, used 
to derive the travelling wave solutions from the standing waves. This symmetry is absent 
for the CH-NLS equation, as it is for the CH equation. In order to see this fact, one can 
easily compute the time derivative of the centre of mass C := J x\m\‘^dx and obtain 


d 

dt 


C = i 



uux)dx. 


(5) 


The right hand side is not a conserved quantity, thus the evolution of the centre of mass 
does not follow a simple linear equation in the time variable, as for the NLS equation, 
where the right hand side would be 2P. 


On the non-integrability of the CH-NLS equation. First of all, the CH-NLS equation 
is only known to admit a weak ZCR, as defined in [1]. This equivalent formulation of 
Equation Q cannot guarantee its complete integrability. It is actually not yet known 
if the CH-NLS equation is completely integrable. We briefly expose here the argument 
against its complete integrability. The possible second Hamiltonian structure which 
could give a bi-Hamiltonian formulation for the CH-NLS equation would be 


0 ~ 
*(1 - 0 


( 6 ) 
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Together with P and M, they will produce the travelling wave equation rut = m^. Then, 
the well-known theorem of Magri would apply to produce the hierarchy of integrable 
equations and there associated conserved quantities. However, this argument fails 
because the two Hamiltonian structures J and K are not compatible. This failure arises 
in the Jacobi identity oi J + K and does not seem to be easy to £x by modification of the 
Hamiltonian structure^ We can therefore not predict the existence of other conserved 
quantities, such as the Hamiltonian itself. 


3. Standing wave solntions 


We derive here a one parameter family of localised standing wave solutions of the CH- 
NLS Equation ([^. Although this solution cannot be promoted to a travelling wave 
solution because of the lack of a Galilean transformation, the standing wave solutions 
still exhibit interesting properties such as peaked profiles. This type of standing waves 
has already been discovered and investigated by [51 [H] in the context of short pulses in 
optical hbres. 

We recall the usual form for standing waves where (j){t),a{x) are real functions to 
be determined 

u = (7) 


thus 


m = {a{x) — and Ut = icftu. 

Inserting this ansatz into Equation Q and using the separation of variables yield 

~^^xx (n Oi axx')(^^ ^ n,j,) I . 


a — a^ttxx 

This allows us to set (p = ct and integrate the resulting equation after its multiplication 
by Ox to obtain 

+ (1 + 2 a^(c — a^))a^ — a^( 2 c — a^) = 0 . 

We set the integration constant to 0 as we will only look for solutions vanishing at 
±oo. The travelling wave ODE can then be derived after calculating the roots of this 
polynomial for Ox and reads 
1 




■(1 J- 2a^(c — a^)) ± 1^(1 J- 2a'^c)‘^ — 


( 8 ) 


The height of the stationary wave is given as a function of c by Omax = found by 
setting Oa: = 0 in Equation (|^. This calculation neglected the fact that there is a square 
root in the equation. Taking it into account gives another limit for the height of the 
solution and its maximum gradient 

1 -h 2a^c 


Q-v 


2a 


®x,max ^^4' 


(9) 


§ We are indebted to A.N.W Hone and J.P. Wang for spotting the non compatibility of the two 
Hamiltonian structures. 










On a deformation of the nonlinear Schrodinger equation 


6 


We thus have two different standing wave heights which corresponds to two types of 
solutions. If Omax < ^max ^^e solutiou will be smooth, otherwise a peak will appear at 
the maximum. Indeed, the condition ^ = 0 is equivalent to having a jump in Ox, thus 
a peaked profile. The critical value for c for which the peak appears is Ccut = and if 
a —?• 0 then Ccrit —t oo, which is compatible with the non-peaked NLS standing waves. 
We display on Figure 1(a) the phase portrait of the standing wave solution, given by 




(a) Phase portraits of standing waves. 


(b) Standing waves 


Figure 1. The left figure displays the phase portrait for different values of c for 
a = 0.5. In this case, the critical value of c for the transition between smooth and 
peaked profiles is Ccrit = 2, plotted in magenta. The right hgure displays the standing 
wave profile for the same values of c. For large values of c, the peak becomes sharper as 
the jump in the derivative of a becomes larger. As compared to a standard hyperbolic 
secant profile plotted in dashed line, this standing wave solutions always have sharper 
peaks and larger tails. 


the right hand side of Equation ([^ and on Figure 1(b) the standing wave solution for 
the same values of c. We fixed a = 0.5 and varied c such that the two different solutions 
appear. These solutions have also been tested in the numerics and remain stable, even 
with the peaked profiles which require higher spacial resolutions. 


4. Solitary wave solutions 

The main class of travelling wave solutions that will interest us here are the solitary 
waves. In the case of NLS, the solitons have the form 

■u(x,t = 0) = e*‘''^y4sech(y4a;) (10) 

which can then be transformed under a Galilean boost to give the full time dependant 
solution. As noted before, this Galilean boost does not exist for the GH-NLS equation 
which makes this method useless here. An important point regarding this solution is 
that A and v are two independent free parameters which define a two parameters family 
of solutions. In the GH-NLS equation, the Helmholtz operator will always couple the 
phase and the amplitude of the solution such that the shape of the solitary wave will 
depend on the speed. The definition of A and v must therefore be adapted if one want 
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to characterise the CH-NLS travelling wave solutions. This fact makes the calculation of 
the travelling waves difficult with the usual travelling wave ansatz and we will therefore 
only focus here, as a first step, on numerical solutions. We want to observe solitary 
waves of the CH-NLS, hence we did some numerical experiments with the hyperbolic 
secant as initial data. 



Figure 2. Numerical solutions of the CH-NLS equation are displayed in both panels. 
The left panel displays a selection of stable solitary waves. They are not travelling wave 
as they periodically oscillate in time. The right panel displays the space-time plot of 
a larger initial condition which undergoes fission and produces two solitary waves as 
well as dispersive radiations. 


Recall hrst that the CH-NLS equation contains a length scale a, which can be 
converted into a scale in wavelength given by 1/a. Provided that the spectral width 
of the initial data is small compared to 1/a, namely that |ko ± Afi:| 1/a, where Ak 
is the spectral width of the solution and kq its maximum wavelength, the solution will 
behave as a NLS soliton. In contrast, the higher order terms of the CH-NLS equation 
will start to play an important role and other type of solutions and behaviour will be 
observed, such as breathers type solutions. Note that in the case of our initial condition 
in Equation (10) we have the simple relations Ak = 1/A and Kq = c. The numerical 
experiments displayed on Figure 2(a) are snapshots at f = 150 with initial profile given 
by Equation (10) with parameters A = 1 and u = 1,2, 3,4, 5. We also used a = 0.3 such 
that the numerical scheme could converge with realistic resolutions. In each case, one or 
more solitary wave will emerge from the initial profile as well as dispersive radiations on 
both sides. In order to illustrate this effect, we plotted on Figure [2(b) a simulation from 
a larger initial pulse with A = 1.5 which undergoes an initial fission into two solitary 
waves and dispersive radiations. In most of our simulations, the solitary waves at large 
times are not pure travelling waves, but are of breather type, namely they exhibit 


periodic motions in time. They contain substructures easily identihable on Figure 2(a) 
which periodically oscillate as the solitary wave moves. Another interesting fact is 
that the four solitary waves displayed in this figure have approximately the same phase 
gradient (roughly corresponding to a i/ in Equation ([l0|)), which makes the breathing 
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phenomenon as well as the speed of the solitary wave dependent on the shape of the 
envelop of the solitary wave and not just on the spectral peak position, as for NLS. 
Indeed, the sharper the solitary wave, the faster. Although the solutions of Equation 
Q are very rich, we will not further pursue this investigation here as we are missing 
the analytical solutions of the equation and because precise numerically based results 
about the properties of the solutions are out of the scope of this short study. 


5. Collision of solitary waves 


Another interesting phenomenon of soliton equations that will be explored in this section 
is the interaction of solitary waves. For integrable equations it is guaranteed that the 
solitons will keep their identity after the collision but is therefore not sure for the 
CH-NLS equation. Nevertheless we numerically observed collisions with this property, 
despite some numerical errors (see 


Appendix A 


(a) Slow and large solitary wave, (b) Slow and narrow solitary (c) Fast and large solitary wave. 

wave. 

Figure 3. These figures display three different collisions between two solitary waves 
of the CH-NLS equation, extracted from initial value problems with sech profiles. The 
left solitary waves have a narrower spectrum than the other two, hence have a collision 
more similar to the NLS collision. The other two collisions have higher amplitudes 
and sharper structures than the first, as one can see from inspection of the differences 
in colours between the solitary wave and the structures of the collision. This is due 
to higher wavelengths in the solitary wave spectrum, thus higher nonlinear dispersive 
effects. The numerical scheme is even not able to fully capture the interaction and 
high frequency waves remain at the location of the collision. We refer to [Appendix A| 
for more details on this numerical error. 



We have run three different collisions, with solitary wave extracted from initial 
value problems as in the previous section. We used three different types of solitary wave 
for this experiment. The first solitary waves in Figure 3(a) [ have a narrow spectrum near 
the zero wavelength. The second solitary waves in Figure 3(b) are sharper than the first, 
thus have a larger spectrum. The last solitary wave in Figure 3(c) have a similar profile 
than the hrst, but with a higher speed. The hrst collision is therefore similar to a typical 
NLS soliton collision whereas the others differ. The spectrum of last two type of solitary 
waves contain higher wavelengths, closer to the scale 1/a, thus the effect of nonlinear 
dispersion is more apparent than for the first collision. The main difference is that the 
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amplitude and sharpness of the structures created during the collision are bigger than 
for a standard NLS collision. This corresponds to the excitation of higher wavelengths 
than for NLS equation despite the regularisation feature present in the linear dispersion 
relation, as mentioned in the introduction. This shows that the nonlinear dispersion of 
the CH-NLS is an important factor for such collision behaviours, mainly used to describe 
rogue wave events. 


6. Modulational instability 


We will end this study with some considerations on the modulational instability (MI) 
usiin], fundamental mechanism of the focussing NLS equation and which is also present 
in the focussing CH-NLS equation. We will follow the standard derivation of the MI 
gain for the NLS equation. We first perturb the plane wave solution of Equation ([^ 
with e <C /3 as follow 

u{x) = u^{x,t) = {(3 + e(a;, 

and then obtain the linearised equation for e 


iet — iof^^xxt + ^xx + 2/3^(e 4- e) — 0. 


( 11 ) 


The growth rate of the MI denoted g can be found using e of the form e = e+e 




+ 


e_e 


iojt—iKX 


and is given by 


g{n) = Real 


1 -f 




( 12 ) 


The regularisation comes only from the convolution terms of the linear dispersion, hence 
the MI behaviour will be similar to the NLS equation, but with a smaller growth rate. 
This effect is typical of nonlocal interactions in NLS equations, see for instance 


We plotted on Figure 4(a) the gain g for different values of a with /3 = 1. The effect 
of the nonlocality only reduces the gain, not the size of the instability region as in [TH] . 
After that the instabilities have developed, nonlinearities take over and the solution 
evolves toward a train of pulses. It happens that rogue wave events sometimes emerge, 
with the Peregrine solution being one of the most frequent. We refer to [12] for a 
recent numerical study of the appearance of rogue waves in the NLS equation. We have 
performed similar experiments with the CH-NLS and extracted candidates for Peregrine 
solutions at different values of a. All the simulations we did started with a plane wave 
of amplitude 1 and of two periods in a periodic domain of length 50. The instabilities 
developed naturally from the noise induced by the numerical errors. We observe a 


delay in the instability growth, as predicted by Equation (12) and the formation of 


rogue events. We displayed the best candidates of Peregrine solution in Figure 4(b) for 
a = 0,0.12,0.18. Similarly to the simulations of collisions of solution, the nonlinear 
effects are stronger than for the NLS and numerically difficult to capture, hence the low 
value of a that we could use. The shape of the solution are similar and close to the 
Peregrine solution, but the amplitudes are higher with higher values of a. The standard 
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(a) Modulational instability gain. 



(b) Peregrine solutions. 


Figure 4. The left panel displays the modulational instability gain of Equation (121, 
derived from the linearisation of the CH-NLS around the plane wave solution with 
P = 1. The gain is reduced but the unstable region not. The right panel exposes a 
selection of Peregrine type solutions of the CH-NLS equation, extracted from numerical 
simulation of modulational instabilities with the same colours as for the left panel. The 
nonlinear dispersion of the CH-NLS increased the size of this extreme event, making 
it interesting rogue wave candidates. 


factor of 3 for the ratio of the amplitude of the Peregrine solution with the background 
amplitude is obviously recovered with a = 0 and reaches values as large as 6 with only 
a = 0.18. 


7. Summary and open problems 

We presented here our latest results on the CH-NLS Equation ([^, a deformation of the 
NLS equation, previously derived in [1]. We hrst exposed the Hamiltonian structure, 
the conserved quantities such as the mass and momentum of this equation and discussed 
its lack of complete integrability. Further investigations of the integrability of the CH- 
NLS are needed, but would have to face technical difficulties, owing to the Helmholtz 
operator. 

We then studied a class of solutions which do not move, the so called standing wave 
solutions. We explicitly found the phase space ODE but only studied the standing waves 
in the case where the wave vanishes at inhnity, thus the non-periodic case. The main 
Ending was the existence of a transition form smooth to peaked solutions by varying the 
only parameter of this class of solutions. Unfortunately the lack of Galilean invariance 
of the CH-NLS equation does not allow the standing wave solutions to be promoted to 
travelling waves. The analytical solution for the travelling wave is actually difficult to 
compute, hence we performed numerical experiments in order to find and study them. 
From initial conditions given by the NLS solitons, we observed the emergence of solitary 
waves and dispersive radiations. Most of the solitary waves exhibit a periodic motion 
in time, similar to the evolution of a breather of the NLS equation. The structure of 
these oscillations were not studied in details as this would require analytical solutions 
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or more systematic numerical studies. 

We then studied the collision of solitary waves with the help of numerics and found 
that the general behaviour is similar to the NLS collisions but with sharper structures 
of higher amplitudes. This fact is interesting as it is linked to the study of rogue 
waves, extreme event which are known to appear in nature and are usually described 
by equations of NLS type. Apart from collisions of solitary waves, rogue waves can be 
found by studying the nonlinear regime of the modulational instability, which is still 
present in the CH-NLS equation. We also found rogue wave events of Peregrine type, 
but sharper and with higher amplitudes than NLS Peregrine solutions. 

We hnally want to mention that in most of our numerical experiments, the solutions 
converged well when increasing the space and time resolutions, except at the location of 
the extreme events such as the collisions of solitary waves or the Peregrine solutions. We 
did a brief analysis of this fact, but we let the development of more accurate numerical 
schemes for future works. 
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Appendix A. Numerical scheme 


The numerical integration of Equation ([^ was performed with a pseudo-spectral 
numerical scheme. We used a Fourier spectral method in space, with second order 
centred hnite difference approximation for the space derivatives in the nonlinear term. 
The time integration is second order implicit for the linear part and first order explicit 
for the nonlinear term. The general scheme is given by 


Ui+l 


Ui 


dt 


1 -|- 


2 '^i+i “1“ I r n 
-K ---+ {m(| 


u\ 




)}; 


(A.l) 


where u{k) is the Fourier transform of u{x) and i runs through the timesteps. Although 
the numerical scheme for the linear term only is not dissipative, the nonlinear term 
will introduce dissipation. For each simulation we ensured that the spatial resolution 
is good enough the capture the full spectrum of the solution. We noticed that the 
CH-NLS equation with large a is numerically difficult to solve, owing to the nonlinear 
dispersion. In the simulations of collisions and modulational instabilities we observed 
the appearance of sharper peaks which eventually breaks into a noisy high frequency 
set of waves. These waves do not move, as predicted from the linear dispersion relation, 
and are purely dispersive. In Figure ^ we plotted the errors on the conservation of 
mass and momentum during the collision displayed in Figure 3(b) for three different 
timesteps. This illustrates the first order of convergence of our scheme as well as the 
difficulties to properly resolve the interactions of solitary waves. 
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(a) Numerical error of mass M. (b) Numerical error of the momentum P. 

Figure Al. Both panels of this figure show the errors for the collision displayed in 
Figure |3(b)[ The error for the mass is relative whereas the error for the momentum 
is absolute, because the value of the momentum is close to 0. As expected, the 
convergence is of order 1, but the collision is never completely resolved, as seen from 
the jumps in the error of both the momentum and mass. 




























